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In Paper I in this series we constructed evolution equations for the complete gauge- invariant linear 
perturbations of a time-dependent spherically symmetric perfect fluid spacetime. A key application 
of this formalism is the interior of a collapsing star. Here we derive boundary conditions at the 
surface of the star, matching the interior perturbations to the well-known perturbations of the 
vacuum Schwarzschild spacetime outside the star. 



I. INTRODUCTION AND OVERVIEW 

In Paper I in this series Q| we constructed evolution equations for the complete gauge-invariant linear perturbations 
of a time-dependent spherically symmetric perfect fluid spacetime, with arbitrary two-parameter equation of state 
p = p(p, s). We isolated true degrees of freedom for all perturbations. They obey a hyperbolic system of wave and 
transport equations. A key application of this formalism is the interior of a collapsing star. Here we derive boundary 
conditions at the surface of the star, matching the interior perturbations to the well-known perturbations of the 
vacuum Schwarzschild spacetime outside the star. Our formalism is now complete and ready for numerical work. 

As in Paper I, we use the covariant and gauge-invariant perturbation formalism of Gerlach and Sengupta and 
our notation is compatible with theirs. The combination of Paper I and the present Paper II is intended to be self- 
contained, but we refer the reader back to Paper I for some definitions and results. In our presentation, we go from 
the general to the specific. 

In Section [n| we discuss continuity conditions across a hypersurface in spacetime where the stress-energy tensor 
is finite but possibly discontinuous. We argue that the appropriate choice of continuous perturbation fields are the 
perturbations of the 3- metric and extrinsic curvature of the surface with contravariant indices. We give those quantities 
in terms of metric perturbations and their first derivatives. 



In Section [II we then restrict the background spacetime to be spherically symmetric, but still allow for arbitrary 
matter, and for the matching surface to be cither timclike or spacelike. We identify a complete set of continuous gauge- 
invariant perturbations in the notation of Paper I, separately for the axial and polar perturbations, decomposing all 
tensor quantities into components in a frame aligned with the surface. This simplifies and corrects results of P,H. 



In section [V we restrict consideration to perfect fluid matter and a timclike matching surface, namely the surface 
of the collapsing star, and state the matching conditions across the stellar surface in terms of variables specialized to 
fluid matter. The conditions for the axial perturbations are fairly simple. 

In section [v] we bring the continuity conditions for the polar perturbations into a final form that shows clearly how 
information crosses the stellar surface in both directions ("extraction" and "injection"), in a way that is natural for 
a numerical evolution. 

We now sketch the form our main results are going to take. In Paper I we characterized the axial perturbations in 
the interior by a tangential fluid velocity perturbation (3 that obeys an autonomous transport equation, plus a metric 
perturbation n that obeys a wave equation with (3 as a source, n is defined both in the interior and exterior of the 
star, and obeys the same equation in the interior and exterior. (In the exterior the matter source (3 simply vanishes.) 
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P requires no matching because the stellar surface is formed by fluid worldlines and is therefore a characteristic of the 
transport equation for (3. The matching condition for II is simply that it is continuous at the surface. 

In Paper I we also characterized the I > 2 polar perturbations by fields x, k and -0 that admit free Cauchy data 
and evolve autonomously. By definition, these are all metric perturbations, but we showed that while \ characterizes 
gravitational waves, k characterizes the sound waves, and ip characterizes a second tangential fluid velocity perturba- 
tion. (The matter perturbations properly speaking can be reconstructed from these primary variables by quadratures 
if they are needed.) Consequently, \ an d k obey wave equations and ip a transport equation (all of which are coupled). 
It is well known that the / > 2 polar perturbations in the exterior are gravitational waves characterized by a single 
field Z that obeys a wave equation (the Zerilli equation ||). ip is transported along the fluid, and therefore, like the 
perturbation variable /3, requires no explicit matching. The variables x, k in the interior, and Z in the exterior, how- 
ever, obey wave equations. This means that these three variables can and must be specified on a timelikc boundary. 
The matching problem for the polar perturbations therefore consists in finding \ an d k just inside the surface from 
Z and its derivatives just outside (injection), and finding Z just outside from x, k, ip and their derivatives just inside 
(extraction) . 

Matching the polar perturbations in the special cases I = and 1 = 1, which do not admit (completely) gauge- 
invariant perturbations and therefore require (partial) gauge-fixing, is also discussed. 

In the appendixes we compare our results to previous work, and give some intermediate steps of the calculation. 



II. JUNCTION CONDITIONS ON A HYPERSURFACE 



A. Choice of continuous perturbation objects 

Given a spacetime (M 4 , g^v) containing matter, we want to derive junction conditions on a hypersurface S where 
the stress-energy tensor is allowed to be discontinuous. Our principal application will be the surface of a star, with 
fluid matter in the interior, and vacuum in the exterior, but for now we still allow the matching surface to be either 
spacelike or timelike. We assume that the metric is at least twice differentiable on each side of the boundary, but we 
do not assume that one coordinate patch covers both sides of the boundary. Therefore we do not assume that the 
metric components or their derivatives are continuous. 

The condition that the stress-energy does not have a <5-function singularity on the surface (a "surface layer") 
translates into the condition that the induced metric and extrinsic curvature of the surface are the same on both sides 
of the surface ||. These tensors are 

iptu = g^v T n^n v , (1) 

^■^iv — TlfAiai v — "F ^/ij'CK^ ^i/j (2) 

where n M is the unit vector orthogonal to the surface, with n^n^ = ±1. Here and throughout this paper the upper 
(lower) sign applies when is spacelike (timelike). (Therefore the upper sign will apply when we restrict to a stellar 
surface later.) The semicolon denotes a covariant derivative with respect to the metric g^ v . The tensors and 
are symmetric, and are intrinsic to the hypersurface in the sense that their contraction with n M vanishes. 

It is helpful to think of the two sides of the surface as two distinct spacetimes, each with a boundary. The two 
halves match without a surface layer if the 3-tensors and are equal. This in turn means that one can locally 
find 3 coordinates on each boundary so that the components of and are equal. To ask, for example, if is 
continuous is not meaningful. 

Clearly, the matching conditions for linear perturbations of the spacetime will be the continuity of the linear 
perturbations of e^ v and , but there is a complication: because the metric itself is perturbed, the perturbations 
of a tensor with its indices in different positions are no longer the same. For the metric itself, for example, with the 
notation 

A(g^) = K»> ( 3 ) 

and the convention that all indices are moved with the background metric g^ v and its inverse g^ v , one has 

A(.gH = -hf". (4) 

Similarly, for the perturbation of any tensor t^ v , 

A(V) = 9u. P 9u\A(t px ) + h^u + t^hxu- (5) 



2 



Which of A(e M „) and A(e A " y ), and of A(i pu ) and A(i pv ), should be made continuous? We are not aware of an argument 
in the literature, and therefore give one here. 

Let the surface be defined by the level surface / = of a scalar field /. The unit normal vector is then the 
normalized version of the gradient 

<^(±/, A /,^)- 1/2 , n M = n" = «Tn„. (6) 

The tangent vector to any curve in the surface / = gives us a vector field X p that is by definition intrinsic to 
the surface. The gradient field f t/l is normal to any of these vectors in the sense that X^f^ = 0, and this fact is 
independent of the metric. Similarly, a contravariant tensor is intrinsic, independently of the metric, if and only 
if / )M = on any of its indices. By extension we may define i M „... = g^ a gup ■ ■ • i a/3 " to be intrinsic, but this 
definition of an intrinsic covariant tensor depends on the metric. Therefore, it is clear that the perturbations A(z A " y ) 
and A(e A " y ) must be continuous (and so A(i M1/ ) and A(e MJ ,) are not). 

In defining linear perturbations of a spacetime by subtracting fields at a point in the perturbed spacetime from the 
corresponding fields at a point in the unperturbed spacetime, the well-known problem arises that there is no unique or 
preferred map that identifies points in the unperturbed and perturbed spacetime. Changing the map (infinitesimally) 
changes the (linear) perturbations, even if the gauge has been fixed in the background. The point identification map 
is in practice provided by introducing coordinates on the perturbed and unperturbed spacetimes, and the linear gauge 
freedom in the perturbations then arises as the absence of a unique or preferred coordinate system on the perturbed 
spacetime, even if the background has a preferred coordinate system adapted to its symmetries. The gauge freedom 
in the linear perturbations arises over and above the gauge freedom in the background. 

The surface of a star has a coordinate-independent significance, and it is therefore natural to identify the perturbed 
with the unperturbed surface. In fact, subtracting a background spacetime from the perturbed spacetime makes sense 
only if the point identification map maps the perturbed surface to the unperturbed surfaces: otherwise we would have 
to subtract, for example, a point in the fluid interior of the perturbed star from a point in the vacuum exterior of the 
background star. The resulting density "perturbation" would then be the background density, which in general is not 
small on the surface. 

To avoid this problem, we choose a perturbation gauge in which the perturbed surface coincides with the background 
surface, obtain the matching conditions, and then transform back to the perturbation gauge in which we actually 
want to work. To formalize this, let us call AX the perturbation of a tensor X in a completely arbitrary gauge. By 
Lie-dragging AX along the vector field 

e ee ±^A/n^. (7) 

we obtain the same perturbation in a new, auxiliary, gauge: 

AX ee AX - C^X (8) 

In this gauge the surface has not moved, because we see that A/ = 0. We call this the surface gauge. (This gauge 
is of course not unique. Rather out of the four gauge degrees of freedom in a completely arbitrary gauge, the three 
parallel to the surface remain free.) In the following, we carry out some intermediate steps in surface gauge, and 
indicate this by using an overbar. As long as our final result contains only gauge-invariant fields, it does not matter 
in which (partial) gauge we have obtained it. 



B. Continuous perturbation objects in terms of metric perturbations 

In surface gauge, with A/ = 0, the perturbation of the normal vector n M is 

A(n M ) = ±-n x n p h Xp n ll , (9) 

so that 

A(n") = ±^n x n p h Xp n" - h pv n v . (10) 

Note that A(n„) is proportional to n^, while A(n p ) is not in general proportional to n p . This is another reflection 
of the fact that the natural intrinsic objects are contravariant: if their contraction with vanishes, so does their 
contraction with n M + A(n ([t ). For the perturbation of a covariant derivative we use the formula 
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A{X^, V ) = [A(X»)]. U + A(T»xu)X x (11) 
and its obvious extensions to other tensors, where the "perturbation of the Christoffel symbol" is a tensor: 



Ar'V = -g»P (h p x-u + h pv ,x - h Xv ;p) ■ (12) 



The continuous perturbations in surface gauge are 

A(^) = -i^ a i vfs h a p, (13) 
A(eH 



h a ,3 - i^i^-ihxo;? + h X p; a - KP;X). (14) 



We can replace continuity of A(e'- iV ) by continuity of the shorter expression 

A(e^) - e</A(0 - e a »A{i» a ) = ±^n a n^h aP - i^ a i^n x ^(hx a -,i3 + W - Kp-,x)- (15) 

A(i^„) and A(e /J1/ ) are not continuous, but one finds by explicit calculation that 

i^A(i a0 ), i» a i^A{e aP ) (16) 

are continuous. (Gerlach and Sengupta ||,|| take the continuity of these expressions, or rather their equivalent in a 
general gauge, as their starting point.) 

III. PERTURBATIONS OF SPHERICAL SYMMETRY 
A. Spherical background 

We now restrict ourselves to a spherically symmetric background, but still allow the matching surface to be either 
timelike (upper sign in all following equations) or spacelike (lower sign). As in Paper I, we write the spherical 
background metric in a 2+2 covariant decomposition as 

g MV = diag (gAB,r 2 -f a b) , (17) 

where j a b is the unit metric on the two-sphere, and the background stress energy tensor as 

= diag (t AB , Qr 2 j ab ) . (18) 

Let n A be the vector field normal to the matching surface, with length squared ±1. Let u A be the vector field 
tangential to the surface, with length squared =pl. The 1+1 metric gAB can be written in terms of this orthonormal 
basis as 

9AB = ±{-u A u B + n A n B ). (19) 

We also use the notation g A B\c = for the covariant derivative in two dimensions, v A = r~ 1 r tAj and the frame unit 
derivatives / = u A f^ A and /' = n A f tA . 

The induced metric and extrinsic curvature of the matching surface of the background spacetime are 

ipv = diag(+u j4 u B , r 2 7 ah ), (20) 
e^y = di&gi^v u A u B , Wr 2 j ab ). (21) 

Therefore the scalars r, v = n A \ A and W = n A v A — r' jr are continuous. Note that although v contains a derivative 
of n A , it depends only on the surface itself. As r is continuous everywhere on the surface, its unit derivative f along 
the surface, and therefore U = u A v A , is continuous. From the continuity of f and r' follows the continuity of r \ A r ,A 
and hence of the Hawking mass m. From the Einstein equations we have that n A n B t AB and u A n B t AB are continuous. 
The list of continuous quantities can be extended by taking dot-derivatives of continuous quantities. The quantities 
that are not required to be continuous include u A u B t ABl Q, the Gauss curvature 1Z of g AB , and /i = u A \ A - Note that 
at the surface of a perfect fluid star (where u A coincides with the fluid 4-velocity) u A n B t AB and Q vanish identically. 
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B. Axial perturbations 



In the following we use the notation of jjj for the gauge-dependent and the gauge-invariant matter and metric 
perturbations, without repeating all the definitions. That notation uses the general metric and stress-energy pertur- 
bations of Gerlach and Sengupta B, adding only specific notation for perfect fluid matter. Our notation here is also 
compatible with the notation of Gerlach and Sengupta in their papers on the matching conditions Note that 

for the axial perturbations we do not need to use surface gauge, as A/ is a polar perturbation. 

The general junction conditions for axial perturbations in the framework of Gerlach and Sengupta are continuity of 

I > 1 : u A hf la \ (22) 
I > 2 : h, (23) 



from continuity of the induced metric (|13|) , and 



I > 1 : u A n B (h a £H - + 2v A h% xiai ) , (24) 



l A\B U B\A 

l>2: n A (h a / lal - h, A ) , (25) 



from continuity of the extrinsic curvature (|15|) . 

In the axial sector we have four gauge-dependent continuity conditions (for I > 2) and one gauge freedom. Therefore 
we find three gauge-independent continuous objects: 

1>1: tt = e AB (r- 2 k A )\ B , (26) 
I > 2 : n A k Al u A k A . (27) 

II is of particular interest as it obeys a wave equation with source terms given purely by matter perturbations. The 
Einstein equations allow us to reconstruct k A later from II and the matter perturbation L A , as 

Z>1: {l-l){l + 2)k A = 16nr 2 L A -e AB (r 4 n) lB . (28) 

Here e AB is the totally antisymmetric covariant tensor with respect to g AB . It can be given in terms of the basis as 

e AB = ±{n A u B - u A n B ) (29) 

We can use this to translate the continuity conditions on k A into conditions on II, the only dynamical variable we 
work with. As II is continuous at all times at the matching surface, its derivative along the unit vector u A in the 
surface, EI, must also be continuous. Multiplying ( p8| ) by n A and using ( |27j ) when I > 2, we find continuity of 

Z > 1 : n A L A . (30) 

This condition must be obeyed automatically if we have chosen the matching surface consistently with the matter 
equations of motion. Multiplying ( p8[ ) by u A we find that 

I > 1 : IT - lQnr~ 2 u A L A (31) 

is continuous, where we have used continuity of r'. This is a condition that needs to be imposed on the initial data 
for II and the matter. It is compatible with the matter conservation equation. The continuity of II itself needs to be 
imposed during the time evolution as a dynamical boundary condition. 

For I = 1 the variable II is still gauge-invariant, but it no longer obeys a wave equation. Instead it is determined 
by the matter perturbations through 

1 = 1: r 4 II = 16ttT, where r 2 L A = e AB T |B . (32) 

The continuity of II implies the continuity of T, and (|3^) and ([H]) still hold. For I = 0, there are no gauge-invariant 
axial perturbations. Our results coincide exactly with those of Gerlach and Sengupta in [jj). 
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C. Polar perturbations 



The general junction conditions for polar perturbations in surface gauge are continuity of 

I > : ii = u A u B h AB , (33) 

~i 2 = K, (34) 

1>1: J 3 = u A h A oU \ (35) 

I > 2 : z 4 = G. (36) 

from continuity of the induced metric, and 

Z>0: ei=u A u B n c (h CA \ B + h CB \ A -h AB \ c ) +n C \ c n A n B h AB , (37) 

e 2 = n A - K, A - 2v A K - l SLtll h P/^ T n c v c n A n B h AB , (38) 

/ > 1 : e 3 = ^ s (S^ - /^ ar + - 2v A h B °^) , (39) 

/ > 2 : e 4 = ti a fe olar - \r 2 G >A - r 2 v A Q\ (40) 



from continuity of the extrinsic curvature. Note that these are all scalars, and from their continuity at all time follows 
that of their derivative along the surface, e.g. of u A K >A = K. 

Now we go back to the general gauge in which A/ is arbitrary. The continuous fields i\ to e 4 in general gauge 
are i\ to e 4 plus terms proportional to A/ and its derivatives. Next we find those linear combinations that are 
gauge-invariant for I > 2. As a rule of thumb, this is done by eliminating G and both components of the vector 
p A = h A olal — (l/2)r 2 G\ A . From this argument we would expect "8 — 3 = 5" continuous gauge-invariant fields. But 
we note that 

N = <pAf -p A n A = <p(Af -p A l A ) , (41) 

is the gauge-invariant perturbation of the scalar / times the background quantity (p. Therefore N is gauge-invariant. 
In surface gauge it reduces to —n p A — — e 4 , which is continuous. Therefore we have 6 rather than 5 continuous 
gauge-invariant perturbations. (JV is the same variable as the N of Gerlach and Sengupta Q.) 

N describes the deformation of the surface: while it is formally gauge-invariant in the bulk, it has a physical 
meaning only on the surface. To understand its significance better, we note that the perturbed / vanishes on the 
perturbed surface, so that 

[/ + A/] (x» +Ax») = Af + f ifl Ax» + 0(A 2 ) = 0, (42) 

so that, in any gauge, 

N = -tpAx A f, A -p A n A = -n A (Ax A +p A ). (43) 

The "surface displacement" Ax M is a gauge-dependent perturbation that transforms as vector field. Regge- Wheeler 
(RW) gauge is by definition the gauge where p A — 0, so we can characterize — N as the normal proper displacement 
of the surface in RW gauge. Note that this is a gauge-dependent statement. In comoving gauge, for example, the 



surface is by definition not displaced. 
The six continuous fields are 

I > : h= u A u B k AB + 2vN, (44) 

I 2 = kT 2WN, (45) 

E x = T2TZN + vn A n B k AB + u A u B n c k AB \ c - 2 (u A n B k AB - nY - 2v 2 N, (46) 

E 2 = ±Wn A n B k AB -k'T 8mL A u B t AB N ± (u 2 - 3W 2 + N T 2U (u A n B k AB - N^j , (47) 

I > 1 : E 3 = u A n B k AB -2N + 2UN, (48) 

l>2: E A = N. (49) 



These fields are continuous in any gauge and for the values of I indicated. They are gauge- invariant for I > 2, but 
only partially gauge-invariant for 1 = 1, and not gauge- invariant for / = 0. 
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IV. PERFECT FLUID MATTER 



We now specialize to perfect fluid matter, and to the case where n A is spacelike. We therefore take the upper sign 
in the equations above. In the stress-energy tensor, we have 

u A u B t AB = p, n A n B t AB = p, u A n B t AB = 0, Q = p. (50) 

The surface on the star is defined by p — 0, and p is therefore continuous, but p can be discontinuous. The fluid 
four-velocity is tangent to the surface, and therefore we use the notation u A for both the unit tangent vector to the 
surface, and for the fluid four-velocity inside, and we use n A both for the vector that is normal to the surface, and 
normal to the fluid four-velocity in the interior of the star. As in paper I, we choose n A to point outside and u A to 
the future. / must then increase with radius to obtain /' > in accordance with the definition (^J). 

A. Axial perturbations 

In the interior, the dynamical degrees of freedom are a matter velocity perturbation [3 and a metric perturbation 
II. In order to work with fields that are O(l) at the origin, we rescale them as 

(3 = r l+1 p, U = r l - 2 Il, (51) 

(3 obeys an autonomous transport equation (it is, of course, transported along with the background fluid), while n 
obeys a wave equation with (3 as a source term (Eqs. (65) and (71) of Paper I). In the exterior, the matter perturbation 
(3 is not defined, but the metric perturbation II is, and it obeys the same field equation just without its (3 source term. 
We can therefore work with the same equations of motion in both the interior and exterior, with [3 defined to vanish 
identically in the exterior. (3 parameterizes tangential fluid motion, and therefore there is no reason why it should 
vanish just inside the surface. It will therefore be discontinuous at the stellar surface. This discontinuity is consistent 
with the equation of motion, as (3 is transported along fluid worldlines, and so could be discontinuous between any 
two fluid world lines. 

Continuity of IT and II is equivalent to continuity of H and n, as r and f are continuous. For perfect fluid matter, 
the stress-energy perturbation L A is given by 

L A = (3(p + p)u A , (52) 

and so L A n A vanishes identically. It is therefore automatically continuous, as suggested above. Continuity of (|3l]) is 
equivalent to continuity of 

II' + 16nrpf3 (53) 

in terms of the rescaled variables, where we have used continuity of II and r' . 

Matching for the axial perturbations is straightforward. The continuity of (|53| ) is a constraint on the initial data S, 

n and (3. It is conserved by the evolution equations. (Note that II' is generally discontinuous if p ^ on the surface.) 
During the evolution, one has to impose the continuity of H. The case I = 1 requires no matching at all, as II is not 
a dynamical variable. 

B. Polar perturbations 

The variable N is not defined in the bulk but only on the boundary. As we want matching conditions for the bulk 
variables, we must eliminate N from the matching conditions. Here we do this for perfect fluid matter. 

The stellar surface is defined by p — 0, and so it is natural to identify / with — p, with the minus sign chosen so 
that / increases with radius. It does not matter that p is not defined in the exterior, as we shall only need to calculate 
N just inside the surface. So N is related to the pressure perturbation, which in turn is determined by the density 
perturbation u> and entropy perturbation a through the equation of state. In appendix |^, we derive the following 
result: 

(54) 

just inside the surface 



> o 



N 



Ca 
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Note that if c 2 vanishes on the surface (as it does for example for polytropic equations of state) ui diverges while N 
is finite. 

A second approach to calculating N is to relate the time derivative of the radial surface displacement TV to the 
radial fluid velocity perturbation 7. In appendix we show that this leads to the expression 



I > 



N-fiN 



7+ 2 



(55) 



just inside the surface 



We can use this equation as a better behaved alternative to ( p4|) by integrating it as an ODE in time along the surface. 

If one combines the dot-derivative of (|54|) with (55), one does not obtain an identity, but a new boundary condition. 
Using the equation for &, and the fact that c 2 and C are constant on the surface, it can be written as 



/ > 



fi + - j (c 2 s lu + Co) + (7 + ^ ) {v + Cs') on the surface. 



(56) 



This equation relates u> to 7, while the matter equation of motion in the bulk relates u> to 7' (Eq. (99) of Paper I). 
Using the bulk equation for tb, we obtain a relation between 7 and 7'. This new equation can be used as a dynamical 
boundary condition governing the reflection of sound waves at the boundary. But its origin is the purely kinematic 
boundary condition that the Lagrangian pressure perturbation Ap vanishes at the surface. 

We can now address the continuity conditions for the bulk variables. We use the Einstein equations to replace TZ. 
We also express the tensor Uab in terms of its components in the fluid frame. Finally, we assume that the stellar 
surface is timelike, and so pick the upper sign. The continuous perturbation fields then become 



I > 



I > 1 
I > 2 



h 
h 

E 2 

E 3 
£4 



X + k - 2r] + 2uN, 
k - 2WN, 

4m 

8n(p + p) 3- 



N + v( X + k) + ( X + k - 2ry)' + 2(iip + 2^ + 2N - 2v 2 N : 
1±Z(Z + 1) 



W(x + k)-k' - 8ir P N + [ U 2 - iW 2 + 

-ip - 2N + 2UN, 
N, 



N + 2U 



U + n) 



(57) 
(58) 

(59) 

(60) 

(61) 
(62) 



where N is given in terms of the matter perturbations by either ( |54|) or (55). Note again that while these fields are 
continuous in any gauge, for the values of I indicated, they are completely gauge-invariant only for I > 2. We now 
consider the three cases I = 0, I = 1 and I > 2 separately. 



1. I > 2 



For perfect fluid matter, rj = for Z > 2 by virtue of one of the perturbed Einstein equations. Taking linear 
combinations of the six continuous fields I\ to £4 and their dot-derivatives, we obtain a simplified set of continuous 
gauge-invariant fields: 



I > 2 : Ci=N, (63) 

C 2 - k, (64) 

C3 = X, (65) 

C 4 = ip, (66) 

C s =k' + 8npN, (67) 

C 6 = x 1 + 2^. (68) 



Note that the coefficients p and p, can be discontinuous on the stellar surface. Note also that we do not need the value 
of N in the exterior, as its coefficient 8irp vanishes there. 
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2. 1 = 1 



For I = 1, r\ does not vanish, and E4 is no longer continuous. As our variables arc now only partially gauge- 
invariant, we make the partial gauge choice k = 0. In this partial gauge, the following five fields arc continuous, and 
gauge-invariant under the remaining gauge freedom: 

l = l,k = 0: G\=N, (69) 

C 2 = X- 2*7, (70) 

C 3 = ip, (71) 

Ci = W X - 8npN, (72) 

C 5 = (u + W) X +(x- 2r?)' + 2^- (73) 

Note that \ IS n °t continuous for I = 1. 



3. / = 



For Z = 0, both E3 and £4 are no longer continuous. Our variables are now not gauge-invariant at all. We begin 
by fixing the gauge partially, again by setting k = 0. The four continuous (but gauge-dependent) fields in this partial 
gauge can be written as 

I = 0, k = : Ci = N, (74) 

C* 2 = X - 2r?, (75) 

C3 = VKx - 871-piV + 2Uip, (76) 

C 4 - (1/ + VK) X + (x - 27?)' + 2(p + U)il> + 2i>. (77) 

Note that \ an d i> are not continuous for I = 0. 



V. MATCHING CONDITIONS FOR THE POLAR PERTURBATIONS 



The physical I > 2 polar perturbations comprise fluid convection, sound waves and gravitational waves in the 
interior, characterized by variables tp, k and X respectively, but all coupled together. In the exterior there are only 
gravitational waves, and these can be characterized by a single variable Z first found by Zerilli. The matching is 
non-trivial because at the surface Z does not simply coincide with \. 

To achieve clarity, we formulate the matching conditions as two sets of boundary conditions, x, ip and k are 
considered as evolving on a spacetime with a timelike boundary, on which one can freely specify certain boundary 
conditions. The most natural ones are Dirichlet boundary conditions for \ and k, the two variables that obey 
wave equations. Similarly, Z obeys a wave equation on another spacetime with timelike boundary, and one can 
freely specify a Dirichlet boundary condition for it. In this view the matching problem consists in finding Dirichlet 
boundary conditions for Z given Xi ^ and k and their first derivatives in the interior (extraction), and Dirichlet 
boundary conditions for x and k given Z and its first derivatives in the exterior (injection). 



A. Vacuum exterior 



The exterior spacetime in spherical symmetry must be the Schwarzschild spacetime, with metric 

ds 2 = -a- 2 dt 2 + a 2 dr 2 + r 2 dQ 2 , (78) 

where a 2 is shorthand for (1 — 2m/r) -1 . Here r and t are the Schwarzschild coordinates, and m is the mass of the 
star, r is identical with the area radius r, and m is identical with the Hawking mass m, which is constant in the 
exterior. The matching conditions in the background spacetime are continuity of r and m at the stellar surface. As 
in Paper I, we use a hat to distinguish the radial frame from the fluid frame. On the Schwarzschild background, in 
the radial frame, we have 
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U = p, = 0, W=(ar)~ 1 , v = {aT l ) , r = mr^ 2 a. 



(79) 



The frame derivative ' of all these quantities also vanishes. The frame derivatives expressed in Schwarzschild coordinates 

are /' = a~ 1 f. r and / = af it . It is useful to introduce the tortoise radius r* by dr* = a 2 (r)dr, so that /' = a/.r„. 

Zerilli Q first found that the I > 2 polar perturbations of Schwarzschild obey a single wave equation for a single 
variable, with all other perturbation variables obtained by quadratures. A gauge-invariant single variable obeying the 
same wave equation (the Zerilli equation) was constructed by Moncrief J?J . In appendix ^ we rederive the Moncrief 
variable and the Zerilli equation in our framework. The final result is that the variable 



with 



Z = A(r)x + B{r)k + C{r)k' (80) 



2r _ r[a 2 {l 2 + l + -2ar 2 

~ a 2 (Z 2 + / + l)-3' B ~ a 2 (l 2 + l + l)-3 ' a 2 {l 2 + 1 + 1) - 3 ' (81) 



obeys the wave equation 



with the potential 



Z l \A - V{r)Z = -Z+ Z" + vZ' - V(r)Z = 0, (82) 



7( r )- !(fiD 6m r 2 A(A + 2) + 3m(r - to) 
r 2 r 3 (rX + 3m) 2 

where X = (I + 2)(l — l)/2. Note that the Zerilli equation can also be written as —Z M + Z.r^r, ~ a~ 2 V(r)Z — using 

the tortoise coordinate. The initial data Z, Z can be set freely on a Cauchy surface, and Z evolves autonomously. 
The other metric perturbations are given algebraically in terms of derivatives of Z, and obey the Einstein equations 
automatically. 

For I = and I = 1 all exterior perturbations are pure gauge. We review this in appendix H 



B. I > 2 matching 

In the interior the principal part of the free evolution equations is —\ + x" + • ■ ■ VA k + c 2 fc" + . . . ip' and ip. The 
evolution equations therefore have 5 characteristics, namely the light cone, the sound cone, and the fluid 4-velocity. 
Two of these characteristics travel from the outside to the inside, and so two quantities can be determined freely on 
the stellar surface. A convenient choice of these is \ an d k. Similarly, the principal part of the free evolution equation 
in the exterior is —Z + Z", so that we can specify one quantity freely on the stellar surface. It is usefully chosen to 
be Z . The matching problem now consists in determining Z just outside from k and ip and their derivatives just 
inside (extraction), and x an d k just inside from Z and its derivatives just outside (injection). 



In section [V, we gave the continuous quantities C\ to Cq in tensor components in the fluid frame. In order to write 
down equalities between quantities just inside and just outside the stellar surface, we must express these quantities 
in terms of Z and its derivatives. For numerical work one would probably want to use comoving coordinates in the 
star and Bondi or polar-radial coordinates outside. Therefore we use fluid frame derivatives inside, but radial frame 
derivatives outside. 

Frame tensor components and related quantities in any two frames are related by hyperbolic rotations through an 
angle £ (for vectors) or 2£ (for 2-tensors). We use the formulas given in Appendix E of Paper I, and take hatted 
quantities to refer to the radial frame and unhatted quantities to the fluid frame. From (|7^), and Eq. (E4) of Paper 
I, we find that 

sinh£ = aril, cosh£ = arW. (84) 

These expressions are evaluated just inside the surface. 

We write down some matching conditions that we shall need as intermediate results. The continuous fields C\ — N 
and C'2 — k are frame-independent, and therefore they are the same just inside and just outside the surface. The 
tensor components 4> = x + k and ip are also continuous, and we just need to transform then from the fluid to the 
radial frame. The continuity of C3 and C4 therefore gives 
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4> — cosh 2^ — sinh2£ ip, 
ijj = cosh 2^ ijj — sinh2£ (f>, 



(85) 
(86) 



where the left-hand side is just inside, and the right-hand side just outside, k' and k transform by a rotation through 
the angle £. k is continuous because k is, but k' just outside is equal to k' + 8irpN just inside. Putting this together, 
the continuity of C2 and C5 gives 



k' + 8npN = coshf k! + sinh£ k, 
k = cosh£ k + sinh£ k' . 

We obtain the extraction equation from the definition (|80| ) of Z, the inverse of 
It is 



(87) 



361) and the inverse of 



Z = A (cosh 2^ (j) + sinh2£ ip) + (B — A)k + C cosh£ (k' + 8npN) - sinh£ k 



(89) 



where the left-hand side is evaluated just outside and the right-hand side just inside, and where A, B and C are the 
coefficients defined in Eq. (|Sl|). Substituting these coefficients, the expressions (|84|), and the definition <fi = x + k, we 
obtain the final version of our extraction equation 



Z = rk 



2r 4 



(I + 2) (I - l)r + 6m 



(W 2 + U 2 ){ X + k) + 2UWtp - W{k' + 8npN) + Uk 



(90) 



The injection equation for k is obtained from the constraint (All) and continuity of k. The injection equation 
for x is obt ained from the c ontinuity equation ([sBI), the definition (f> = x + k, the injection equation for k, and the 
constraints (A12) and (A13). The two injection equations are therefore 



(All) 
Z 

r 



1(1 + 1) 



1 - 



2m 



Z.r- 



6 m 



(I + 2)(l - l)r + 6m r 



cosh2^ JA12| ) - sinh2^ ( gT| ) + (cosh 2^ - 1) ( gig ) 
2m 
r 



1 



[(W 2 + U 2 ) ( [All ) - 2C/V^ ( gl| ) + 2f7 2 JHl|) l , 



(91) 



(92) 



where the equation numbers in brackets stand for the right-hand sides of those equations, evaluated just outside the 
surface, and where the left-hand sides are just inside the surface. We have not written out ( |92] ) in full because no 
important simplifications occur when one expands it. Equations ( |90f]92| ) are the main results of this paper. 
We have not used the continuity of Cq. 



C. I = 1 



For 1 = 1, the polar perturbation variables we use are no longer fully gauge- invariant. It is necessary to make a 
partial gauge choice, and in Paper I we chose k = 0. This still leaves a residual gauge freedom worth one free function 
of time only. In paper I we fixed this remaining freedom by demanding that at the center the leading order in r of 
the variable r\ vanishes, or fj — 0(r 2 ) in the notation of Paper I. In the exterior we also make the partial gauge choice 
k = 0, and again this leaves a residual gauge freedom parameterized by one free function of time. As reviewed in 
Appendix^, this freedom can be used to set all perturbation variables in the exterior equal to zero, which shows that 
the exterior perturbations are pure gauge. This gauge is not the same as the one fixed by imposing fj = 0(r 2 ) at the 
center. We therefore drop the latter - it was imposed only in the absence of a better choice and does not simplify the 
field equations. The equations of Paper I are valid for any way of fixing the residual gauge freedom, and in particular 
for the gauge choice 77 = at the surface that we adopt now. 

With the complete gauge choice k = everywhere and 77 = at the surface we find that x — "0 = V = in the 
exterior. Transforming this to to the fluid frame, we have x = — V = just outside the surface. From the continuity 
of C2, C3 and C4, we have 
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I = 1, k = 0, 7/ = just outside : ip — 0, (93) 

X = 2r), (94) 

V=^N, (95) 

just inside the surface. These equations can be used as boundary conditions for the integration of the ODE system 
(A12-A14) of Paper I. From the continuity of WC 5 - UC 2 , we have 

r\v\ 2 D( X -2r 1 ) + 8np(v + W)N = 0, (96) 

just inside the surface. (See Paper I for the definition of the radial derivative D.) If we now use Eqs. (A13) and 
(A12) of paper I to eliminate D-q and D\, ( p3[]95| ) to eliminate Xj ip and f], and ( |54| ) to eliminate N, and use p = 
at the surface, we obtain an identity. 

For I — and I = 1 we use the matter perturbations u) and 7 as dynamical variables. Together they describe sound 
waves, and therefore we need one boundary condition at the surface. It is given by (pF 



D. I = 



The situation for I = is similar to the one for I — 1. We have already made the partial gauge choice k = 0. In 
order to fix the gauge further, we impose the further gauge choice ip = —n A u A kAB = 0, where \jj is the equivalent 
of the frame component ip in the radial frame, instead of the fluid frame. (The reason for this choice is discussed in 
paper I.) This gauge choice is equivalent to 

2UW 

^ = 0,k = 0: 4,= __(„_ X ). (97) 

We shall use this to eliminate ip from C3 and C4. 

The gauge choice k — ip — still leaves a small residual gauge freedom. In Paper I we used this last gauge freedom 
to fix fj = 0(r 2 ) at the center, but as for I = 1 this is not the best choice to make in the presence of a vacuum exterior, 
where all perturbations are pure gauge. Instead we fix the gauge so that all perturbations vanish in the exterior (see 
Appendix [§]), so that k = r) = 'ip = x = everywhere outside, and therefore also ip = x = . Continuity of C2 and 
C3 then gives rise to the boundary conditions 

1 = 0,k = 0,i> = 0,1] = just outside : x = 2*7, (98) 

W* + U* 4*p N m 
? > w 2 -U 2 w { ' 

The continuity of C4 gives again rise to an identity, similar to the I = 1 case. Again, we need the matter boundary 
condition (j56|). 
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APPENDIX A: DERIVATION OF THE ZERILLI EQUATION 

Here we only consider the case I > 2. Recall that in the interior, \, ip and k evolve autonomously and can be 
considered as the true degrees of freedom. Three constraints give the matter perturbations a, 7 and lo directly in 
terms of derivatives of x, "0 an d k. In vacuum, however, these matter perturbations vanish, so that the same equations 
become three constraints on the Cauchy data for x, ^> and k. They are 
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= (k)' - Wjc 







1(1 + 1) 



= ij/ + 2vip 4 
The evolution equations in vacuum are 



X - 



X" = -2 



2v 2 - 



l{l + l) 
2r 2 

(X + 

f 2 k. 

6 m 



k) 



(i-m 



2r 2 



W X ' - 2Wk\ 



k = 



-Wx'-Pk'-^ix 



(x + k)- (5*> - 2W)x > + 
(l-l)(l + 2 



(/-!)(! + 2). 



k) 



2r 2 



il} = 2v(x + k)+x' 



(Al) 
(A2) 
(A3) 

(A4) 
(A5) 
(A6) 



Note that in obtaining the vacuum evolution equations from the evolution equations inside fluid matter, we h ave set 
c 2 = 0. Setting c 2 to a different formal value would correspond to adding that value times the co nstr aint (A2) to the 
evolution equation (A5). Here we have chosen to write the evolution equation for k in the form (A5) which does not 
contain k" . 

We note that in vacuum ip already plays a passive role, while x already almost obeys an autonomous wave equation, 
containing only k as an unwanted term. This suggests looking for a variable of the form (80) that obeys a wave equation 
of the form (|82| ) above, to be constructed from equations (A2), (A4) and (A5). Note that the new variable Z should 

not contain k in order to avoid bringing tjj back into the game. 

Introducing the ansatz (^) in equation ( |S2| ) and using the above equations to eliminate second derivatives, we 
obtain 



- Z + Z" + vZ 1 - Vox! - V lX - V 2 k - V 3 k' = 0, 
where the Vi are certain combinations of A, B, C and a. It is clear that Vq should vanish: 



V 



2(aC + ra 2 A)' 



0, 



G 



-raA 



Cl 

a 



(A7) 



(A8) 



where c\ is a constant. The coefficients Vi/A, V2/B and V3/C have to be equal to each other and are all equal 
to the "scattering potential" V . Using our result for C this gives two coupled nonlinear differential equations of 
second order for A and B. With the choice c\ — one of them is first order and can be easily integrated (define 
A„ = — n + a 2 (I 2 + I + 1) for integer n): 



*-%=' [ -A 1 A + 2B? = 0, 



B = ±A 
2 



c 2 , 



(A9) 



A C raA 

where c 2 is another constant. Having expressed B and C in terms of A, the other equation gives a single differential 
equation for A, which again can be easily integrated if c 2 = 0: 



B 



Vi 
.4 



1(1 + 1) a 
Ai A 



A, 



A 



0, 



.4 



rc 3 
A3' 



(A10) 



We fix the overall constant C3 in the definition of Z as C3 = 2, which makes our Z equal to the variable found by 
Moncrief. Using these expressions for A, B, C we finally obtain the coefficients given above in (|8l|). 

The initial data Z, Z can be set freely on a Cauchy surface, and Z evolves autonomously. The metric perturbations 
fc, x an d 4> can be reconstructed from Z using the vacuum perturbative Einstein equations flAl[ ), ( |A2j ) and (A3) given 
in the appendix. When we introduce the expression for Z into these equations the derivative terms cancel out and fc, 
X and if; are obtained as algebraic expressions in Z and its derivatives: 

(All) 



(A12) 




(A13) 
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Note that a' — —a m/r . The remaining perturbative vacuum Einstein equations ( A4), flAq) and (|A6|) are now linear 



combinations of the Zerilli equation and its derivatives. Z is therefore the true degree of freedom in the exterior. 

APPENDIX B: L = 0, 1 EXTERIOR POLAR PERTURBATIONS 

a. I = 1: 

Together with k and ip, now we also have the variable r\. These variables are still invariant under gauge 
transformations generated by a vector field £ A on M 2 , but not under gauge changes generated by r 2 £Y' a with £ a 
scalar field. We impose the condition k = to eliminate this freedom. The equations of motion are then greatly 
simplified: 

X' = --X, k = ~l (Bl) 
r r 

V / = I^!v,, j = ?—l { x- 3ri); (B2) 
ar ar 

??? = _X+(-3 + 2a^ (B3) 



They can be integrated to give 



;(t) ^ = rgtffft) ^ = 2mg(t) - r 3 d t 2 g(t) (B4) 



(r — 2m) 2 ' (r — 2m) 2 ' 6m(r — 2m) 2 

where the free function of time g(t) is a residual gauge freedom, not eliminated by the condition k = 0. We are 
still allowed to perform changes of the form £(t,r) = a 2 (r)f(t), where f(t) is an arbitrary function of time. Under 
that change we find Ag(t) = 6mf(t). We can use this freedom to choose g(t) — 0, so that all our perturbations 
vanish. This can be interpreted as a small displacement of the perturbed star so that the center of mass position is 
not perturbed, and therefore the vacuum exterior is not affected by the perturbation. 

b. I = 0: 

In this case we use the gauge k = i/j = 0. The equations are even simpler: 

X' = -~X, 1 = 0, r/ = 0. (B5) 
r 

Their solution is 

2Am , . 

X = — 5-, V = 9(t), (B6) 
r — Ira 

where the constant Am is a perturbation of the mass of the star and g(t) is an arbitrary function of time which again 
represents residual gauge freedom, now related to time reparameterization. To first order in perturbations, we can 
write the perturbed metric as 

2m + 2Am\ ,. . , . 2 / 2m + 2Am\ 1 , 9 . ,„ , „ , . , 

1 {[1 + g(t)]dt} +1 dr 2 + r 2 (d9 2 + sin 2 6d(/) 2 ). (B7) 



The mass perturbation Am is constant in time. In physical terms this is so because the star can lose or gain mass- 
energy only to second order in perturbation theory. We consider only Am = here, and treat any change in the mass 
of the star as a change in the background solution. 
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APPENDIX C: CALCULATION OF N IN TERMS OF FLUID VARIABLES 



We have derived Eq. (|54|) in two ways. 

1. The surface of a perfect fluid star is defined by the vanishing of the pressure p. We identify / inside the star 
with the negative fluid pressure — p (the sign is chosen so that /' > 0), so that N becomes 

N = -(p' A p, A y 1/2 (Ap - P A P,a) = (p')~ 1 p{c 2 s uj + Ccr ) J ust insidc thc surface. (CI) 

(Do not confuse the pressure gradient p iA with the metric perturbation p A . Note also that p' A p, A = p' 2 because p = 
on the surface at all times, and that p' < 0, which explains the overall sign in the second equality above.) From Eq. 
(44) of paper I we find that 

p' = -v{p + p) (C2) 

and so we obtain Eq. (|54|) above. This derivation holds for all I. 

2. We can also obtain (E3) without identifying / explicitly with —p. The evolution equation for fc, Eq. (88) of 
Paper I, holds on both sides of the surface if one formally sets p = in the exterior. It is therefore continuous. Let 
us first assume I > 2. We bring known continuous quantities such as fc to one side, and find that 

c 2 s {—k" + other metric perturbations) — vh! + 8npCa = continuous. (C3) 

Now we use Eq. (94) of paper I to replace the entire term in round brackets by Sirpu, and we use the continuity of 
C5 to replace fc' by —8irpN plus a continuous term. (It is this step that brings in N.) We obtain 

8ttp(c 2 s uj + Ccr + vN) = continuous (C4) 

The left-hand side vanishes identically in the exterior because p = 0, so by continuity it must vanish just inside the 
surface as well, and we obtain Eq. (p4|). For 1 = 1, Eqs. (88) and (94) of Paper I still holds. We now have r] ^ and, 
by gauge choice, fc = 0, but the final result is the same. In the case I = 0, we can use Eqs. (A18) and (A19) instead 
of Eq. (88) to obtain once again the same result. 

We have derived equation (55) for N in two different ways, too. 

1. The velocity perturbation is the Lie-derivative of the position perturbation along fluid world lines: 

Au" = C u Ax" (C5) 

and so we have, in any gauge, that 

n A A(u A ) = n A {u B Ax A \ B - Ax B u A \ B ) = u B {n A Ax A )\ B - u B n A]B Ax A + Ax B n A]B u A = {n A Ax A J - p,(n A Ax A ), 

(C6) 

where we have used Eq. (42) of Paper I to eliminate n A \ B . A simple calculation using Eqs. (52), (55), (21) and (42) 
of Paper I gives 

n A A{u A ) = 7 + I - (n AP A J + p(n A p A ). (C7) 

Eliminating n A Au A between the last two results, and using (^3|) gives us (|55|). This derivation holds for all /. 

2. Alternatively, we can use the field equations directly. Assume / > 2 for now. Using equation (93) of Paper I, 
extracting the continuous terms and using the continuity of C5 and C5 , we find that 

8?r(p + p) - pN + 7 + -0 = continuous. (C8) 

We use the same argument as above: this vanishes in the exterior, so it must vanish just inside the surface, which 
gives us (pq). For I = 1 and I = the derivation is similar, taking into account rj ^= and fc = 0. 
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APPENDIX D: PERTURBATIONS OF A TIME-DEPENDENT STAR: COMPARISON WITH SEIDEL 



In || Seidel considers axial perturbations of a general time-dependent spherically symmetric star. His variable U 
is —/3 in our notation. Our variable II is related to the interior and exterior variables of Seidel as 

II (our notation) = — — — - (Seidel interior) = — — -—r (Seidel exterior). (Dl) 

v ' 1(1 + l)R 2 1(1 + l)r 3 v ' y ' 

The continuity conditions of Seidel agree with the continuity of II that we find. 

In H Seidel derives evolution and matching equations for the / = 2 gauge-invariant polar perturbations. He uses 
comoving coordinates in the exterior and Schwarzschild coordinates in the exterior. His exterior perturbation variables 
Qi and ip are related to our variable Z, restricted to I = 2, by 

~ 3?77 ~ / 4-7T 

Q 1 = (l + ±-)4Z, ^=^J—AZ. (D2) 

ijj obeys the Zerilli equation. His variable k\ is equal to our k. We have already compared his interior variables with 
ours in Appendix C of Paper I. 

Seidel notes that his interior variable is gauge-invariant when restricted to the boundary and in any gauge where 
the perturbed boundary is at the same coordinate location as the unperturbed boundary (what we call surface gauge) . 
He derives the matching conditions in that gauge. We have 

A~ l q 4 = 2n A p A , (D3) 

with Seidel's notation on the left and ours on the right. This is of course equal to — 2N in the surface gauge A/ = 
in which Seidel works, and is therefore actually gauge- invariant, although Seidel does not make that point. Seidel 
gives an evolution equation for q 4 , Equation (2.89), that can be cast in covariant form in our notation if we assume 
that the energy density (p in our notation, rj = (1 + e)p in the notation of Seidel) reduces to the rest mass density p 
of Seidel when p = 0. The rest mass density p of Seidel is linked to his background metric coefficient A through the 
gauge condition A^ 1 = Anr 2 p. His equation (2.89) is then exactly our equation (|55|). 

Seidel's matching equation (2.84), giving Q\ in terms of interior variables is equivalent to our extraction equation, 
Equation (|90|). His second matching equation (2.86), which gives k in terms of his exterior variable Qi is equivalent 



to our Equation (All). His equation (2.85) for Q\. r is a linear combination of the first two. 



APPENDIX E: PERTURBATIONS OF OPPENHEIMER-SNYDER: COMPARISON WITH CPM 

Cunningham, Price and Moncrief study the axial jl0| (CPM1) and polar |jll[ (CPM2) perturbations of spherically 
symmetric homogeneous dust collapse (Oppenheimer-Snyder collapse). The background solution consists of a spherical 
segment of a dust-filled closed Friedmann solution in its collapsing phase matched to a Schwarzschild exterior. The 
interior metric, in the notation of CPM, is 

ds 2 = -dr 2 + R 2 (t) (dx 2 + sm 2 X dn 2 ) . (El) 

Fluid elements are at constant x- The stellar surface is at x — Xo> corresponding to r — ro(r) = sinxo-R( T )- These 
are just a special case of comoving coordinates, and from Appendix C of [lj we can read off that in this metric, in a 
fluid frame, 

» = U = W »' = u ' = °- ' = H ' = ^T (E2) 

Here the left-hand sides are in our notation. (On the right-hand sides, x is a coordinate, not our perturbation variable 
of the same name.) We also have, from the equation of state p = 0, that 

m = —r 3 p, c 2 s = 0. (E3) 

The variable U of CPM1 is —f3 in our notation. Our variable n is related to the interior and exterior variables of 
CPM1 as 
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IT (our notation) = — (CPM interior) = — -. „ (CPM exterior). (E4) 

1(1 + l)i? 2 sin x l(l + l)r 2 

The continuity conditions of CPM1 agree with the continuity of II that we find. 

We first discuss the polar perturbations of Oppenhcimer-Snyder in our notation. For homogeneous dust, the interior 
equations of motion simplify greatly, to 

- X + X" = ^X + 2W X > + (l + 2) ^~ 1) X, (E5) 
-k = U X + Wk-W X ' + 2(w 2 --^) (x + fc)- {l ~ l ^ 2 + 2) X (E6) 
-'tjj = 2Uip + x . (E7) 

X therefore obeys an autonomous wave equation describing gravitational waves that do not couple to matter pertur- 
bations. The sound wave equation normally obeyed by k has degenerated into an ODE along matter world lines, k 
and ip are therefore obtained by solving ODEs after the autonomous wave equation for X has been solved. 

We now compare with CPM2. On an Oppenheimer-Snyder interior, our perturbation variable x is essentially the 
variable Qi of jO], and is related to their final variable ijj as 

/ 5 „ sin v 

X=^QU X=^, (E8) 

where again the notation on the right-hand side is that of CPM2, where x is a coordinate. We find that our equation 



(E5), restricted to I = 2, is equivalent to the autonomous wave equation (11-41) for the variable tp of CPM2. The 
exterior variables Qi and ■0 of CPM2, are the same as the ones used by Seidel and given above in (|D2|). The exterior 
equation of CPM2 is just the Zerilli equation. 

Although CPM2 can and do ignore k and ip in the interior, they need some information about their values at the 
surface for imposing the matching conditions. CPM2 evolve a variable A along the stellar surface that is related to 



our variable k via 

1 



3 sin x 



(E9) 



Equation (B24) or (II-43c) of |Ll|, an ODE evolution equation for A, is our Eq. (EC), restricted to the boundary. 
Equation (II-43a) of |ll[ agrees with our injection equation for k, Eq. (^lj), if one expresses Z' through Z' and Z . 
CPM2 use it not as an injection equation for k (as already discussed, they do not evolve k at all), but as an extraction 
equation for Z of the form AZ + BZ + CZ' = Dk, using the value of k on the boundary that is obtained from (II-43c). 
The injection equation for X of CPM2 is (II-43b) of pT| . Translated into our notation, is 

3y - 4r ■ 3m m „ 9m k ,„ 

— = Z+ UZ+ ~Z , (E10) 

r 2r + 3m 2r + 3m r 6 2r + 3m r 

where the right-hand side is evaluated just outside the surface, even though the derivatives have been expressed in 
the fluid frame. When we transform the derivatives on the right-hand side into the radial frame, and then used the 

Zerilli equation to eliminate Z, we find that this equation agrees with our injection equation (BS 



APPENDIX F: POLAR PERTURBATIONS OF A STATIC STAR: COMPARISON WITH THORNE ET AL. 

The matching conditions for the I > 2 polar perturbations of Thorne I |Q, equation (19), are corrected in Thorne 
II [^3|, equations (Bla-b), and we compare the corrected version with our results. In the following the left-hand 
sides are Thome's notation and the right-hand sides are ours. We use the fact that our gauge-invariant variables are 
identical with the gauge-dependent ones in Regge- Wheeler gauge, which is the one Thorne et al. use. 

Equation (D4) of our Paper I contains two sign mistakes: the overall signs of 7 + ip/2 and of a should be reversed. 
Correcting this, we see that e-"/ 2 W T , t =r 2 ( 7 +|). O n a static background, this in turn is the frame Wt in our 
notation. Comparing with our equation (^5|), we have Wt = —r 2 N (on a static background only), where we have 
set the integration constant to zero as it does not carry physical information. We also have Hq = H2 = —(x + fc), 
e (\+u)/2jj^ _ ^ an( j g _ _^ Q n tne gtatjc background, e~ A / 2 f >r = f. (Note that Thorne uses /' for / >r , while we 
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use it for the frame derivative along n A .) e x is continuous because the Hawking mass is continuous. Thorne et al. 
also choose e u to be continuous. With these correspondences, we find that the five continuous perturbation fields of 
Thorne II, (Bla-b), are equivalent to the five 

X, k, V, x', k' + 8npN. (Fl) 

These are our continuity conditions, given that fi vanishes for a static background. 

The special case I = 1 has been covered in Thorne V Jl4| . Using the relations (D3-D4) of Paper I we have re-obtained 
all their evolution equations (14) and (15), with the one difference that their (14a) should have an additional term 
—S'Yr/'y (notation of Thorne, see there). We have verified that their regularity conditions (20) at the center agree 
with our conditions (117-120) of Paper I. The matching equations (21b-e) of Thorne V are equivalent to our matching 
conditions (|69|-[73[), if one restricts the latter to a static background and assumes the gauge choice 77 = x — V — 
made by Thorne V in the exterior. 
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